Improved WKB analysis of cosmological perturbations 
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Improved Wentzel-Kramers-Brillouin (WKB)-type approximations are presented in order to study 
cosmological perturbations beyond the lowest order. Our methods are based on functions which 
approximate the true perturbation modes over the complete range of the independent (Langer) 
variable, from sub-horizon to super-horizon scales, and include the region near the turning point. 
We employ both a perturbative Green's function technique and an adiabatic (or "semiclassical" ) 
expansion (for a linear turning point) in order to compute higher order corrections. Improved 
general expressions for the WKB scalar and tensor power spectra are derived for both techniques. 
We test our methods on the benchmark of power-law inflation, which allows comparison with exact 
expressions for the perturbations, and find that the next-to-leading order adiabatic expansion yields 
the amplitude of the power spectra with excellent accuracy, whereas the next-to-leading order with 
the perturbative Green's function method does not improve the leading order result significantly. 
However, in more general cases, either or both methods may be useful. 



PACS numbers: 98.80.Cq, 98.80.-k 
I. INTRODUCTION 

Anisotropies in the cosmic microwave background 
(CMB) radiation and inhomogeneities in the large scale 
structures of the Universe have nowadays become a fun- 
damental tool to study the early universe. This trend 
will persist in the next years due to the data releases of 
Wilkinson Microwave Anisotropy Probe Sloan Dig- 
ital Sky Survey and the launch of the Planck mis- 
sion 0. 

In a primordial inflationary phase, matter and space- 
time fluctuations are generated which inherit character- 
istics from the particular inflationary model studied (for 
an introduction and review see, e.g. Ref. £|). It is there- 
fore very important to be able to determine the power 
spectra of cosmological perturbations for a large variety 
of different inflationary models, so that comparison with 
available and future data will tell us which models sat- 
isfactorily represent the time evolution of the Universe. 
Since exact solutions are not available for cosmological 
perturbations in general inflationary models (except for 
the case of an exponential potential [Bj]) approximation 
methods are very welcome. The method of approxima- 
tion most used is the slow-rollover approximation, intro- 
duced by Stewart and Lyth [(| (see Ref. Q for a review), 
which was inspired by power-law inflation. 

It is also worthwhile testing different approximation 
schemes. In this paper we present Wentzel-Kramers- 
Brillouin (WKB)-type approaches which improve the 



proposals of Refs. jl Ej (see also Refs. El El El El 
Ej, E3 f° r other approximations). For power-law infla- 
tion, the usual WKB method to lowest order reproduces 
exactly the spectral indices of perturbations, but the am- 
plitudes are still rendered poorly. Given the importance 
of the consistency equation between the ratio of tensor to 
scalar amplitudes and the spectral index of tensor modes 
for single Klein-Gordon scalar field inflation 4|, we be- 
lieve that the prediction on the amplitudes should also 
be improved, as is done in Ref. |Ti| by means of the so- 
called "uniform approximation" to next-to-leading order 
(for the details on the formalism, see Ref. [H). 

Let us begin by recalling the Robertson- Walker metric 
in conformal time rj, 



ds 2 = & 2 (?7) [ — d?7 2 + dx ■ dx] 



(1) 



where a is the scale factor of the Universe. Scalar (den- 
sity) and tensor (gravitational wave) fluctuations are de- 
termined on this background by two functions, respec- 
tively denoted as fi = (is ■= aQ (Q is the Mukhanov 
variable E3) an d M = := o,h (h is the amplitude 
of the two polarizations of gravitational waves) 0, ^| , 
which satisfy a one-dimensional Schrodinger-like equa- 
tion 



— + k 2 
drj 2 



M = 0, (2) 



where k is the wave-number. The time-dependent fre- 
quency is given by 
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n 2 (k,ri) = k 2 



(3) 



where z = z§ a<j)' /TL for scalar (0 is the homoge- 
nous inflaton) and z — z^ := a for tensor perturba- 
tions, primes denote derivatives with respect to 77, and 
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H (:= a' /a) is the conformal Hubble parameter. Eq. (j2J 
must be solved together with the condition that the 
modes are initially plane waves for wavelengths much 
shorter than the Hubble radius, 



lim n(k,rf) 



-ik rj 



'2k 



(4) 



where d = d$ '■= 1 and d — c?t ■— n/Stt/topi for scalar 
and tensor perturbations (mpi is the Planck mass). The 
dimensionless power spectra of scalar and tensor fluctu- 
ations are then given by 



fc 3 



Ms 

zs 



4fc 3 



(5a) 



The spectral indices and their runnings are defined at an 
arbitrary pivot scale k„ as 



n s - 1 := 



dlnV c 



d In k 



fe=fc» 



and 



d 2 lnP c 



(dlnfc) 2 



cut := 



dlnfc 



d 2 lnV h 



k=k, 



(dlnfc) 2 



k=k, 



(5b) 



(5c) 



The evolution of the Universe is usually described by 
means of a set of flow equations [H, 1 2(1 l21j . The zero 
horizon flow function is defined by 



H{Nj) 
H(N) 



(6) 



where H := a/ a is the Hubble rate and dots denote 
derivatives with respect to cosmic time dt = a(rj) di]. In 
this expression, N is the number of e-folds, N := ln(a/aj) 
[where a, = a(?7i)] after the arbitrary initial time ti = 
tiirji). The hierarchy of horizon flow functions is then 
defined according to 



tn+l 



din |e,, 
~dN~ 



n > 



(7) 



and we recall that inflation takes place for ej < f . 

The effective potentials U — Us and [/ = Ut can now 
be expressed in terms of e\, e 2 , and only as 



a 2 H 2 

Us(v) 
a 2 H 2 



ei 



(8) 



■ £2 



£2 £3 



and Eq. can be solved exactly for e\ constant in time, 
which is the case of an exponential potential for the in- 
flaton. When the e ra have an arbitrary time dependence, 
approximate methods are necessary. 

A WKB analysis of this problem has been recently 
presented in Ref. to first order in the adiabatic (or 



"semiclassical" ) expansion (see next Section and Ref. |23 
for precise definitions), in which the Langer transforma- 
tion m 



k 

~Ha~ 



£ := In ( ^- | , x:=(l-ei) 1/2 e-*/ 2 ,i , (9) 



introduced in Ref. is used [34| in order to improve 
the accuracy. Such a transformation brings Eq. into 
the form 



d' 

i- it, 

da; 2 



u 2 {x) 



X = , 



(10) 



where the (new) frequency tu(x) in general vanishes at the 
classical "turning point" x = x*. On sub-horizon scales 
(x 3> x*), the mode x oscillates (lu ~ k) and the stan- 
dard WKB approximation can be applied without diffi- 
culty. On super- horizon scales (x <^ x*), the perturba- 
tions exponentially decay or grow and the standard WKB 
approximation can again be applied. One then matches 
the two approximate solutions at x = x* to obtain the 
super- horizon amplitudes at < i, which determine 
the CMB spectra. This procedure opens up the possibil- 
ity of further improving the knowledge of the spectra by 
including subsequent adiabatic orders. 

However, a straightforward generalization to higher 
adiabatic orders is seriously hindered by the lack of a 
general prescription for matching the two WKB branches 
with sufficient accuracy. This problem was thoroughly 
discussed in Refs. |23,|lj|, where it was also suggested to 
replace the standard (plane wave-like) WKB functions 
with Bessel's functions. The latter indeed remain good 
approximations at the turning point, whereas the former 
require the matching with yet another particular solu- 
tion at some (unspecified) points both to the left and 
right of x — x*. We shall follow the proposal of Ref. |2^| 
for a linear turning point and show that the inclusion of 
higher adiabatic orders indeed improve the results, con- 
trarily to what was found in Ref. [25| . We shall also intro- 
duce a new (perturbative) expansion for general turning 
point, which makes use of the Green's function technique, 
and compare the corresponding corrections for the case 
of power law inflation. 

The article is organized as follows: In Section [H] we 
shall review the application to inflationary cosmological 
perturbations of the standard WKB approximation to 
leading adiabatic order || and its shortcomings. In Sec- 
tion II I II we shall describe in full generality a method 
of improving the results to higher orders, both by the 
perturbative expansion discussed in detail in Section llVI 
and the adiabatic expansion reviewed in Section In 
Section IVll we shall apply the improved WKB approxi- 
mation to power-law inflation which, being exactly solv- 
able, will allow us to test the improved WKB predictions 
against exact solutions. Finally, we shall give our con- 
clusions and observations in Section IVHI 



3 



II. STANDARD WKB APPROXIMATION 



where 



In this Section we begin by recalling that the straight- 
forward application of the standard WKB method leads 
to a poor approximation and that the results can be im- 
proved to first adiabatic order by means of the trans- 
formation |J5J, as was shown in Ref. ||- However, the 
extension to higher orders is more subtle, as we shall dis- 
cuss in some detail. 

The WKB approximation is defined by first introduc- 
ing a "small" parameter 6 > in order to perform the 
formal adiabaticfor semiclassical) expansion of the mode 
functions [I^. Consequently, Eq. @ is formally re- 
placed by 



d?? 



(11) 



and, of course, the limit 8 — *■ 1 must be taken at the end 
of the computation. We then denote the leading order 
term in such an expansion of n by 



MwkbO,??) = 



which satisfies the following differential equation 
dij 2 



(12) 



6— + n 2 (k,7 1 )-6Qn(k,7 1 ) 



MWKB = , (13) 



3 (ft' 



4 ft 2 



2ft 



(14) 



We now observe that ^wkb is the exact solution of 
Eq. (jTT|l in the "adiabatic limit" 5 — ► 0, and is expected 
to be a good approximation to the exact /i (in the limit 
S -> 1) if 



A 



ft 2 



< l 



On using the potentials in Eqs. 
H, one finds 



(15) 

in terms of a and 



As.t = 



If) 



sub-horizon scales 



(16) 



super-horizon scales. 



The result for super-horizon scales suggests that a change 
of variable and function are necessary in order to obtain 
better accuracy. 

In fact, on employing the transformation ©, one finds 
(for 5 = 1) the new equations of motion l|ll)|l . which ex- 
plicitly read 



d 2 Xs(z) 


e 2x 


4 




- ei 


dx 2 


_(l-^i) 2 






- ei 






i 




- ei 


dx 2 


(1-ei) 2 


~ 4 


(r 


- ei 



(3-2 £l )e 2 (l-2 ei )£ 2 e 3 (l-4 £l )6 2 
2(l- £l ) 2 2(l- £l )3 4(l- £l )4 

eie 2 e 3 . {2 + e 1 )e 1 e\ 



XS (x) = (17a) 



ei £2 



XtW = . (17b) 



From u>s(x) and w t (ie), respectively given by the expressions in the square brackets of Eqs. (|17af> and (|17b|) . we now 
obtain 



f A s ,T = e- 2 *0(e„)~O 



sub-horizon scales 



A s = — Ui 4 + ei e 2 e 3 + ^ + ) + 0(efj super-horizon scales 



(18) 



A T = 27 (d <3 + ei e 2 e 3 ) + 



super-horizon scales , 



and it is thus obvious that A becomes very small for 
any inflationary potential satisfying the slow-roll condi- 
tions |12| in both limits of interest (but, remarkably, not 
at the turning point). As a consequence, the WKB ap- 
proximation is now also valid on super-horizon scales. 

One might naively expect that the transformation (JSJ 



will improve the WKB approach to all orders (in the 
formal expansion in 5), thus yielding increasingly better 
results. However, since Q diverges at the turning point 
x — x* , whether the function yU-yvKB can still be regarded 
as a solution to Eq. (jl 1|> for x ~ x* in the adiabatic 
limit 5 — s- becomes a subtle issue. It is precisely for 
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this reason that one usually considers particular solutions 
around the turning point (for example, combinations of 
Airy functions for a linear turning point) which must then 
be matched with the asymptotic form (|12|) on both sides 
of the turning point. The problem with this standard 
procedure is that there is no general rule to determine 
the actual matching points (which cannot be x* where Q 
diverges) in such a way that the overall error be small. 
As a consequence, the standard adiabatic expansion for 
5 -C 1 may not lead to a reliable approximation around 
the turning point, which would then result in large errors 
to higher orders and a questionable further extension of 
all the expressions to the interesting case 5 = 1. 

It was in fact shown in Ref. that a simple exten- 
sion to higher orders of the above procedure |2y] actually 
seems to reduce the accuracy with respect to the leading 
order results of Ref. [8j . We shall see in the next Section 
that one can overcome this problem by introducing mode 
functions which remain good approximations also around 
the turning point. 



III. IMPROVED WKB APPROXIMATION 

We shall now describe a method to improve the WKB 
approximation to all orders by following the scheme out- 
lined in Ref. [23| (see also Ref. |13). The basic idea is to 
use approximate expressions which are valid for all values 
of the coordinate x and then expand around them. Two 
different such expansions we shall consider, namely one 
based on the Green's function technique and the other 
on the usual adiabatic expansion |23[ . 

Let us first illustrate the main idea of Ref. j2^| for 
the particular case u> 2 = C [x — x*) n , with C a positive 
costant. All the solutions to Eq. (|1(J|I for x > x* can be 
written as linear combinations of the two functions 



u±(x) 



U)(x) 



J±m [£(x)] 



(19) 



where 



uj{y) dy 



n + 2 



and Jj/ are Bessel functions 28]. Moreover, for a general 
frequency u>, the expressions in Eq. I|19fl satisfy 



dx* 



+ uj 2 (x) — cr(x) 



u± = , 



(21) 



where the quantity (primes from here on will denote 
derivatives with respect to the argument of the given 
function) 



a(x) 



3 (^') 2 



4 uj 2 



2lu 



I \ u 2 



m "47 * 



(22) 



contains the term Q defined in Eq. (|14|) , whose divergent 
behavior at the turning point x — x» we identified as the 
possible cause of failure of the standard WKB approach. 
For a general (finite ) fr equency, which can be expanded 
in powers of x — x* [3|| , 



u> 2 (x) = C (x — x*y 



c « fa 



q>l 



, (23) 



one finds that the extra term in Eq. (12211 precisely "re- 
moves" the divergence in Q at the turning point. In fact, 
the residue 



lim a (x) 



3 (n + 5) c\ 



3c 2 



2 (n + 4) (n + 6) n + 6 



(24) 



is finite |23l |27| . The finiteness of a at the turning point 
is crucial in order to extend the WKB method to higher 
orders (this point was missed in Ref. 25] ) . It is also easy 
to show that, for the cases of interest, \a/ui 2 \ ~ A as 
given in Eq. I|18(l for x — > ±oo, so that the new approx- 
imate solutions remain very accurate for large For 
slowly varying frequencies, with |cr(x)| small everywhere, 
the functions (|19fl are thus expected to be good approx- 
imations to the solutions of Eq. I jlOj l for the whole range 
of x, including the turning point. 

We can now introduce both the adiabatic expansion of 
the previous Section and a new formal expansion of the 
mode functions by replacing Eq. (|1(J|) by 



S — + uj 2 (x) — 6 a(x) 
ax z 



X=-6sa(x)x, (25) 



where 5 and e are "small" positive parameters. We shall 
refer to expressions proportional to 8 n as the n-th adia- 
batic order and to those proportional to e n as the n-th 
perturbative order. 

It is convenient to consider the solutions to Eq. Ij25(l to 
the left and right of the turning point X* separately. We 
shall call: 



(20) 

region I) where uJ 2 > (on the right of with 



(a:) := y/u 2 (x) , 6(x) := / ufi(y)dy , (26a) 



and the solutions to the corresponding homoge- 
neous equation l|21|l are given in Eq. (|19H : 

region II) where to 2 < (on the left of x»), with 

LJn(x) := v/-w 2 (x) , £n(x) := / Wii(y) dy , (26b) 

J X 

and the solutions to the corresponding homoge- 
neous equation (|21|l are obtained from those in 
Eq. (JTSJ by replacing the J v with the I v |28|. 
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Corresponding expressions are obtained for all relevant 
quantities, for example o\ and an, and we shall omit the 
indices I and II whenever it is not ambiguous. 

Although it is possible to expand in both parameters, 
so as to obtain terms of order 6 P e q (with p and q posi- 
tive integers), in the following we shall just consider each 
expansion separately, that is we shall set 5 = 1 in the 
perturbative expansion (see Section llVjl and e = 1 in the 
adiabatic expansion (see Section El . 



IV. PERTURBATIVE EXPANSION 

For 5 = 1, the limit e — > is exactly solvable and yields 
the solutions 119|) . whereas the case of interest (|10|) is 
recovered in the limit e — > 1, which must always be taken 
at the end of the computation. The solutions to Eq. Ij25(l 
can be further cast in integral form as 



X(x)=u(x)-e G(x,y)a(y)x(y)dy 



(27) 



where u(x) = A + u + {x) + A_ w_(x) is a linear combina- 
tion of solutions 119|) to the corresponding homogeneous 
equation l|21jl. and G(x,y) is the Green's function deter- 
mined by 



dx 2 



+ Ld 2 (x) — <j(x) 



G(x,y)=e5(x-y) . (28) 



The solutions of the homogeneous equation (|21|l in the 
two regions are then given by 



ui(x) 



{A+J +m [6(a)] +A_ J-mMx)}} 



: (x) + ui_ (x) 



(29) 



un(x) 



Ux) 



{B+ I +m [fo(a!)] + B- I_ m [£ n (x)]} 



V ^ii (a;) 
= B+ Un+(x) + B_ Un_(x) 



and the Green's functions can be written as 

7T 



Gi(x,y) 



(x - y)ui-{x) ui+(y) 



(30) 



Gn(x,y) 



l 7T 

7!' 



(x-y) un-(x) u u+ (y) 



With the ansatz 
Xi (x) = A + a + (x) u\+ (x) + A- a_ (x) iti_ (x) 

Xu(x) = B + b+{x) u u+ (x) + B- b-{x) un_(x) , 



(31) 



from Eqs. (|29|l and l|30(l . the problem is now reduced to 
the determination of the x-dependent coefficients 

a+{x) := [l + e J+_(x, xi)] + e —— JZ_{x, x r ) 

Aa 



A+ 

i-(x) := [l + e J+_(x*,a;)] + £ J+. 



_i_(x*, xj 



(32) 



&+(x) := [l — iel+_(x, a;*)] — i £ — — J (x, x*) 

-°+ 

6_(x) := [l-iel+_(x f ,x)] - is I+ + (x f , x) , 

in which, for the sake of brevity (and clarity), we have 
introduced a shorthand notation for the following inte- 
grals 



J^(xi,x 2 ) := -j= J aj(y) a w (y) u\ s {y) ui s (y) dy 

= -7= f <*i{v) ^~T\ a ™(y) Jsm i&(y)] J sm Ki(y)] dy 

V3 J X1 uiiy) 

(33) 

If 5 (xi, x 2 ) := -j= J crn(y) & w (y) u ns (y) u n5 (y) dy 

= -7s / a n(y) K(y) ism [61(2/)] hm Kn(y)] dy , 

where the symbols w, s and s are + or — . 

Before tackling this problem, we shall work out rela- 
tions which allow us to determine the constant coeffi- 
cients A± and B± uniquely. 



A. Initial and matching conditions 

We first use the initial conditions i n order to fix the 
constant coefficients A±. 

For the function x such conditions become 



lim x( x ) = d 



1 e l{ x i) -ikrn-Xi/2 



2k 



1- ei(xj) rih(xi ^ Xi/ 2 ^ / 34 s 



where, from the definitions of 77 and x, krji ~ — i£i(xi). 
We then recall the as ymp totic form of the Bessel func- 
tions for x — > +00 [23, |23 |33 1 



J±m [flfa)] 



tt£i(x) 



(ti(x)T^m- J) , (35) 



from which we obtain the asymptotic expression of the 
mode function in region I, 



6 



Xi 



(xt) ~ \J l ^ Xt) e ~ x '/ 2 { e +^> [e- * A + + e"' - (A_ + e A_ J+_(x», a*) + e A+ J+ + (x*, x,)) 

+' (l+f ) + e +' (i-?) 1 (i_ + £ i_ J + -_(i t , + e A + J+ + (i» , ^))] } • (36) 



The initial conditions (|34|l therefore yield 



A, 



-A 



k i- e 2im7r 

6 " ~t~ £ 

1 + eJ+_(x*,Xi) 



(37a) 



(37b) 



Note that A + does not depend on e, therefore the result 
(|37a|l holds to all orders. 

We next impose continuity at the turning point in or- 
der to determine the coefficients B±. For this we shall 
need the (asymptotic) expressions of Bessel functions for 
x — > x* [l^ll^ HU- In particular, 



J±m [6(a)] 



61 (a;) 



1 



r(l±m) 



1 



r(l±m) 



for x — > xj~ 



for a; — > x* 



(38) 



where Y is Euler's gamma function. Near the turn- 
ing point, we just keep the leading term in the expan- 
sion l|23|) . so that 



cj(x) ~ vC I a; — x*|" 



(39) 



£(x) ~ 2 m v C |x — x* I 



We thus obtain the following forms of the mode functions 
near the turning point, 



Xi(x ~ x*) ~ A + a+(x*) 



+A_ a_(x*) 



/2m 



r(l + m) 
^2 m [my/C 



r(i -m) 



-m 

— (40a) 



^2 m 



Xn(x ~ x*) ~ 5+ &+(x*) 



r(i + m) 

/ 2 m [my/C 



r(i -m) 

Continuity across x* then implies that 
£>+ &+(&*) = — A + a+(x*) 

B- &_(x») = A_ a_(x*) . 



-m 

— . (40b) 



(41) 



From the definitions 1M2I) at x = x», Eq. I|37b[) and the 
relations l|41|l we therefore obtain 



B+ _ i | _ c+ 1 "+£j| + ( a: »,x,) 
A+ 1 + e J+_(x*,Xi) 



J — (x^ , xi) 



£ J j (x* , Xj) 



(42) 



B. 



1 + e J+_(x*,Xj) 



1 — is I + _(xf,x*) 



1 -l - £ — (x* , Xi ) 

1 — ie/^_(x/,x») ++ 



(x/,x») 



It is apparent from the above derivation that, since 
the matching between approximate solutions in the two 
regions is always performed at x = x„, no ambiguity 
hinders the evaluation of the coefficients B± and the ac- 
curacy of the approximate solution is therefore expected 
to increase with increasing order. 



B. Recursive (perturbative) relations 

We can finally obtain expressions for the coefficients 
a±(x) and b±(x) by solving the integral relations 13211 . 
We remark that the r.h.s.'s of such equations contain 
a±(x) and b±(x), which can therefore be determined re- 
cursively by expanding them (as well as any function of 
them) in e, that is 

a±(x) = l + ^e 9 <4 9) 

9>1 



(43) 



b±(x) = l + J2e q b± 



(9) 



It will also be useful to define the following integrals 
J™ (9) (xi,x 2 ) := -y= J ai(y)a$(y)ui s (y)uis(y)dy 

(44) 

7™ (9) (xi,x 2 ) := -j= J a u (y) b^ ] (y) u IIs (y) u Hs (y) dy , 

where again w, s and s = ±, and q is a non- negative 
integer. 
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1. Leading order 
We begin by considering the leading (e°) order 



al\x)=bf{x) = l 



(45) 



The corresponding solutions are simply given by 
Xi(x) ~ ui(x) (46a) 
= A 4 



Xn (x) ~ un(a;) 
= -i4 



{j +m [£i(a;)] - e-" m *J_ m [£i(a;)]} 



(46b) 



/ jn(g) 



{/ +m [&(»)] +e +<w " r i'-T n [&(a:)]} 



where we used the zero order forms of Eqs. (|3T 
and gH, 



A_ = B- = -A+e +lm7r 
B, = —Aa 



(47) 



H- > 



and A + is given in Eq. i|37a[l . 

We know that these expressions are also good solutions 
of Eq. near the turning point and, deep in region II, 
we can use the asymptotic form of I v for x ~ Xf — » 

-oo mm, 



-f±m [£ll( 3 



1 



(48) 



On neglecting the non-leading mode e few, we then ob- 
tain 



k 



■Ktxm\ x \aHJ [1 - ei(r/)] o;ii,s(fc,r?) 



,(0) 



n s - 1 ~ 3 + 2 



d£n,s 



d In k 



k=k. 



a s ~2 



d 2 6i,s 



(d In*)" 



(50a) 
(50b) 

(50c) 



k=k. 



and 



T- / /- \ :; e 2fn, T (fc,>7) 



16 iF / jfe 



ttto^ ^ai/y [1 - £1(77)] w nT (A;,r7) 



(0) 



n T ~3 + 2 



d&i; 



d In k 



k=k, 



ax — 2 



d'Cii/ 



(d In jfc)' 



(50d) 



(50e) 



(50f) 



k=k. 



where we have transformed back to the original variables 
k and ?y, and all quantities are evaluated in the super- 
horizon limit (i.e. for k <C aH). 



2. Next-to-leading order 

We now insert the leading order expressions (|45|l into 
the integral relations (|32|l to compute the next-to-leading 
order expressions 

<&\x) = jf_{x, Xi ) - e +4m7r jW(x,Xi) 



a^> (x) = 4°2 (x* , x) - e~ l rn * J™ (re. , z) 
6^(x) - -i lf } _(x,x*) - ie +imn I { °l(x,x*) 



irmr t(0) 



(51) 



6 W (s) = -i J|"i (x f ,x)-ie~ im * r+i (x f , x) , 
where we used the simplified notation J^ ^ = J S g 



„• 7-(°) 



— i m 7r r(0) 



e and I c e 



J s / ^ =: Igg , and also expanded the 



coefficients A- and _B± to zero order in e. 

From Eq. 1)48(1 ■ again neglecting the non-leading mode 
e-fiitx)^ we i-,t a i n the mode function at a; = Xf to first 
order in s, 



~ -A+ V 7 e«"^/) . (49) xn(i/) 

y/2irijJii{xf) 

If we now use Eqs. (I5a|) . (|37a|l and 149|) . and the values 
of d, we recover the results of Ref. |8j for all quantities 
of interest, i.e. power spectra, spectral indices and a- 
runnings 



A+ e^ 11 ^) 

y/2lTLJu(Xf) 



1+e 



+i e 



lfl(x f ,x*) -I^i(x f ,x*) 



-ice 



-ee ' 



J -\ — 7 %i ) "I - 5 



jP(i„a:0 + 4 i(a;.,a:i) }• , (52) 



which, using Eqs. (|5aJ> . 1)3 7a|) . 1(521) and the values of d, 
yields all quantities of interest to first order in e. For the 
spectra we obtain 



1 + 9™^) 



(53a) 



8 



where the relative corrections to the leading order ex- 
pressions (|5Ua|l and (|50d|) are now given by 



(53b) 



+ 



n/3 



I { +i{x, x») - 2 J™ (x, x*) + I™ (x, x*) 



r(0) 



(0) 



S,T 



in which we set e — 1 and n = 1 as required, and write 
S and T to recall the use of the corresponding frequen- 
cies. The expressions for the spectral indices and their 
runnings to this order can finally be derived from the 
definitions (|5b|) and l(5cjl. 



3. Higher orders 



The above procedure can be extended to all orders. 
The coefficients up to and b± must be used to com- 
pute the integrals J™ s ^ and which determine a± +1 ^ 



and b 



(9+1) 



However, general expressions soon become 



very involved, and we shall test the effectiveness of our 
method by applying it to a case of interest in Section lVll 



V. ADIABATIC EXPANSION 

Let us now apply the usual adiabatic expansion with 
the assumption that the leading order be given by the 
functions l|19|) [rather than the more common expres- 
sion l|12|)]. This leads one to replace 



£0*0 



V5 



(54) 



and consider the forms 21 



U± =F±(x; 5) u± (x) + G± (x; S) u'± (x) , (55a) 

where the x-dependent coefficients F and G must now 
be determined. The particular case of a linear turning 
point |29j | (precisely the one which mostly concerns us 
here), can then be fully analyzed as follows. 



A. Recursive (adiabatic) relations 

The x-dependent coefficients F and G can be expanded 
in powers of 5 as 



F(x;<5)=5» (j 



j>0 



(55b) 



j>0 



Upon substituting into Eq. Q25[) with e = 1, one there- 
fore obtains that the coefficients <pu\ and are given 



recursively by the formulae 



f 



70)0*0 



(55c) 



2cj(x)J ( CT(2/) [ 2 ^-i)^) + ^- 1 )^). 



d/y 



where it is understood that the integration must be per- 
formed from Xst to x in region I and from x to x* in 
region II. Let us also remark that, for ui 2 ~ (x — x*) near 
the turning point, the above expressions are finite pj^ . 
whereas for more general cases one expects divergences 
as with the more standard WKB approach. 



1. Leading order 

The lowest order solutions H19fl are correctly recovered 
on setting 



0(0) = 1 and 7(o) = , 



(56) 



in the limit 5 = 1. The relevant leading order perturba- 
tions are therefore obtained on imposing the initial condi- 
tions 1341) and matching conditions between region I and 
region II at the turning point, which therefore yield for 
Xi and xii the same linear combinations U\ of Eq. I|4tia|) 
and tin of Eq. Ij46b|l previously obtained. 



2. Higher orders 

The condition (|56l) makes the formal expressions of the 
coefficients 4>(i) and 7(i) particularly simple, 



0(1)0*0 



7(i) 



2ui(x) J uj(y) 



7('i)(2/) + <r(y) 7(i) (y) 



dy 



(57) 



and the first order solutions are then given by linear com- 
binations of the two functions 

U±{x) = [l + <50 (1) (x)] u±{x)+6 1{1] (x)u' ± (x) . (58) 

It is now important to observe that the initial condi- 
tions l|34|l and matching conditions at the turning point 
to first order in the adiabatic parameter S yield different 
results with respect to the perturbative expansion in e. 
First of all, the coefficient A+, which was left unaffected 
by the expansion in e [see Eq. (|37a|) ]. acquires a correc- 
tion. The relations between A_, B± and A + instead 
remain those given in Eq. (|47|l to all orders, whereas 
in the perturbative expansion they were modified [see 
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Eqs. i|37b[) and 142(1 ]. To summarize, to first order in 6, 
one obtains 

+ ~ V k 1 - e 2im7T 



-A+ e +im7r 



and the perturbation modes deep in Region II (at x <C 
a:*) are finally given by 



A. 



(59) 



B± = TA± , 



Xu(x) ^ ujj(x) <l + S 



0ii(i) {%) ~ 7n(i)(s) ^n(^) 



^ii(g) 
2u>n(x) 



01(1) (^i) + 7i(i)0 E i 



iui(xi 



.(60) 



r 



We note that, in order to obtain the above result, one 
must first take the asymptotic expansion of the functions 
u± and then take the derivative, since the reverse order 
would lead to larger errors. The corrected power spectra 
are then given by 



P c ~7f> [l + g$ a (x f ) 

V h ^vf [l + gf» T (x f ) 
in which we set 8 = 1 at the end, and l3£ 

9(i)s,t( x ) = 2 



(61a) 



0ii(i) 0) - 7ii(i)(» u u (x) + 



, 7r(i)(si) , , \ 
+ ^ 0i(i) 



(61b) 



S,T 



where the indices S and T recall the use of the corre- 
sponding frequencies. As usual, the corrections for the 
spectral indices and a-runnings then follow from their 
definitions (|5b(l and f5c| . 

Given the complicated expressions for u>, it is usually 
impossible to carry out the double integration that yields 
analytically, let alone higher order terms. One must 
therefore rely on numerical computations. In Section lVll 
we shall, in the context of power-law inflation, compare 
the results obtained in next-to-leading order from the two 
expansions so far described in full generality. 



with (3 < —2. The case (3 = —2 is special, since it corre- 
sponds to the de Sitter space-time with constant Hubble 
radius equal to £q . The frequency is given, both for scalar 
and tensor modes, by 



2 „2i 



w\x) = (1+/?)' e 



0- 



(63) 



and the horizon flow functions read 
2 + 



ei 



1+/3 



, n > 1 



(64) 



For the particular case f3 = — 2, we plot the frequency in 
Fig- G] the function £(x) in Fig. [3 and the perturbation 
a(x) in Fig. 




VI. APPLICATION TO POWER-LAW 
INFLATION 

In this model the scale factor is given by 

a(v)=Mv\ 1+P , (62) 



FIG. 1: The frequency lo 2 for (3 = —2. The turning point is 
at x, = ln(3/2) ~ 0.4 



Eq. J2J with the initial conditions Q can be solved 
analytically and the exact power spectra at Xt — > — oo 
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FIG. 2: The function £ for (3 




FIG. 3: The "perturbation" cr for (3 = -2. 



are given by (see, e.g. Refs. |5l l8l l30l l31ip 

1 1 



C 



7r ei m P1 Zq 



(65a) 



7rm P1 «q 



where 



/(/?) := - 

7T 



r (1/3 + 1/21 



2/3+1 



(65b) 



and /(/? = — 2) = 1 40] . The spectral indices and their 
runnings can also be calculated from Eqs. i|5b() and l |5cjl 
and one finds n$ — 1 = tit = 2/3 + 4 and as = ax = 0. 
Exact scale-invariance is obtained for f3 = —2. 



of Section lead to the same result as in Ref. 
1 1 



V h 



7T ei m pl Iq 



16 1 



2 . g (°)(/3)fcW 



(66a) 



7rm pl tg 

where the function g( ' (/3) is given by 

2p 2 ' 3 + 1 
JO)'- ze 



|2/3+l| 



2/3+2 



(66b) 



This yields a relative error for the amplitude of the power 
spectrum 



100 



/(/?)- <? (0) (/3) 



/(/?) 



(67) 



which decreases for increasing \/3\, as can be seen from 
Fig. 01 but is rather large (about 10%) for the de Sitter 
space-time. The spectral indices and their runnings are 
instead predicted exactly as ng — 1 = ut = 2 (3 + 4 and 
Q!g = ax = 0. 



B. Next-to- leading order results 

We shall now compare the corrections to the ampli- 
tude of the power spectra coming from the two different 
expansions we described in Section ITTT1 



Perturbative expansion 



From the next-to-leading expressions in Section flV B 21 
we obtain that the relative correction to the power spec- 
tra is given by the function g™ EEN := Srnf EN = 

■9 (iff™ in Eq ' ^ 53b ^ evaluated at x = Xf < a;*. This 
quantity can be determined numerically and one finds 
that it approaches an asymptotical finite value for Xf — > 
— oo. Some examples are given in Table[I]for Xf = —13, 
from which it is clear that the improvement over the lead- 
ing order is small. Once the correction has been included, 
the total relative error on the power spectra (denoted by 
Ap ) therefore remains of the same order as the leading 



order error A 



(0) 



2. Adiabatic expansion 



A. Leading order results 

The leading order expressions in the perturbative ex- 
pansion of Section ITVl as well as the adiabatic expansion 



We analogously evaluate the coefficients 4>(x) an d 7(i) 
numerically. For the case j3 = —2, the result is plot- 
ted in Fig. |5j from which it appears that the functions 
(f)(1) (x) and 7(i) (a;) converge to finite asymptotic values 
for x — > ±oo. Since the same behavior is also found for 
all (3 < —2, we exhibit the values of such coefficients 
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p 


-2 


-3 


-4 


-5 


-6 


-7 


-8 


-9 


-10 


GREEN 

9(1) 


+9.6- 10" 5 


+7.3- 10" 4 


+7.1 ■ 10 -4 


+6.4- 10" 4 


+5.8- 10" 4 


+5.2 • 10" 4 


+4.8- 10 -4 


+4.4 • 10~ 4 


+4.2 ■ 10~ 4 


A (0) 

iAp 


10.4% 


6.4% 


4.6% 


3.6% 


3.0% 


2.5% 


2.2% 


1.9% 


1.7% 


A (l) 


10.4% 


6.3% 


4.6% 


3.6% 


2.9% 


2.5% 


2.1% 


1.9% 


1.7% 



TABLE I: Next-to-leading order improvement for the power spectra 
The latter is essentially the same as to leading order (Ap'). 



?(if EEN ) and total final error (Ap 1 ') to first order in e. 



p 


-2 


-3 


-4 


-5 


-6 


-7 


-8 


-9 


-10 


01(1) 


+3.2 • 10~ 3 


+1.1 ■ 10~ 3 


+6.0- 10" 4 


+3.6 • 10~ 4 


+2.4- 10~ 4 


+1.7- 10~ 4 


+1.3- 10~ 4 


+1.0- 10~ 4 


+8.1 ■ 10~ 5 


7i(i) 


-2.4- 10" 19 


-5.9 • 10" 20 


-2.6 ■ 10~ 20 


-1.5 ■ 10" 20 


-9.4 ■ 10~ 21 


-6.6- 10" 21 


-4.8 • 10" 21 


-3.7- 10" 21 


-2.9 ■ 10~ 21 


0ii(i) 


+4.5 ■ 10~ 3 


+1.6- 10~ 3 


+8.2- 10~ 4 


+5.0- 10~ 4 


+3.3- 10~ 4 


+2.4 • 10 -4 


+1.8- 10~ 4 


+1.4 • 10~ 4 


+1.1 • 10~ 4 


7ii(i) 


-3.5 ■ 10~ 2 


-1.2 ■ 10" 2 


-6.4- 10" 3 


-3.8- 10" 3 


-2.6 ■ 10 -3 


-1.8 • 10~ 3 


-1.4 ■ 10~ 3 


-1.1 • 10" 3 


-8.6- 10" 4 


9(D 


+1.1 ■ 10" 1 


+6.3- 10~ 2 


+4.5 • 10~ 2 


+3.5- 10~ 2 


+2.9 ■ 10~ 2 


+2.4 ■ 10~ 2 


+2.1 ■ 10~ 2 


+1.8- 10~ 2 


+1.6- 10~ 2 


A (l) 

iAp 


0.83% 


0.50% 


0.35% 


0.27% 


0.22% 


0.18% 


0.16% 


0.14% 


0.12% 



TABLE II: First order coefficients and 7(i), correction (<?(i) ) and total relative error (Ap') for the power spectrum to first 
order in <5. 



10- 



-10 



p 



FIG. 4: The total error on the power spectrum to leading 
order and to next-to-leading order for the perturbative cor- 
rection (superposed upper lines) and to next-to-leading order 
for the adiabatic correction (lower line). 



(with Xi = 20 and Xf = —13) in Table HT1 together with 
the corresponding next-to-leading order relative improve- 
ment [g^?(xf) as defined in Eq. (|61b(l ] and the total er- 
ror on the power spectra (Ap"*). On comparing with the 
last row in Table [IJ we can therefore conclude that the 
adiabatic corrections are significantly better than those 
obtained from the perturbative expansion and, in fact, 
yield the amplitudes with extremely high accuracy (see 
also Fig. In fact, this method seems (at least) as 
accurate as the next-to-leading order in the uniform ap- 



proximation employed in Ref. 



g ore 

m 



VII. CONCLUSIONS 

In this paper we have improved the WKB approxima- 
tion for the purpose of estimating the spectra of cosmo- 



-0.02 - 



-0.03- 



-10 -5 



5 10 15 20 



FIG. 5: The functions (j)^(x) (solid line) and 7(i)(a;) (dotted 
line) for f3 = -2. 



logical perturbations during inflation. We found general 
formulae for amplitudes, spectral indices and a-runnings 
of the fluctuations to next-to-leading order both in the 
adiabatic expansion of Ref. |23] and a new perturbative 
expansion which makes use of the Green's function tech- 
nique. 

We have then applied our method to power-law infla- 
tion in order to test it against exact results. It is known 
that the spectral indices and their runnings are obtained 
exactly in leading order, hence we focussed on the spec- 
tra to next-to-leading order. We have found that the 
perturbative corrections remain too small to yield any 
significant improvement, whereas the adiabatic expan- 
sion to next-to-leading order reproduces the exact am- 
plitudes with great accuracy (see Ref. 0] for a similar 
result with the uniform approximation). This result does 
not however mean that the perturbative expansion will 
not lead to significant corrections in different inflationary 
scenarios. 

One way of understanding the difference between the 
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two expansions in the power-law case may be the follow- 
ing. It is known that the Born approximation in quantum 
mechanics (analogous to our perturbative method with 
the Green's function) is good for high multipoles (angu- 
lar momenta) I of the expansion in spherical harmonics. 
In our approach, the parameter (3 plays the role of the 
multipole index I for the energy levels of the hydrogen 
atom, as can be seen on comparing the frequency (IC33|> 
with the expression given in Ref. [23] (see also Ref. y|). 
The approximation with the Green's function is there- 
fore expected to yield more significant corrections for 
large values of |/3|, which coincide with the regime of 
fast-roll. Indeed, the leading order becomes more and 
more accurate for increasing Since the interesting 
regime for inflation involves small values of (3 ~ —2, it 
is instead the adiabatic approximation which seems bet- 
ter since the horizon flow functions evolve slowly and/or 
the states are quasi-classical 01 ■ This may be the rea- 
son whereby the adiabatic approximation works better 
in the next-to-leading order. However, in more general 



cases, either or both methods may contribute significant 
corrections. In this respect, let us remark that the adi- 
abatic expansion can be straightforwardly applied only 
for a linear turning point |42j, whereas the perturbative 
Green's function method is not so restricted. 

Having assessed the accuracy of the method, it is now 
natural to use it in order to improve current estimates for 
inflationary models which are not exactly solvable. We 
are in fact analyzing the slow roll approximation with 
our method and will report about it in a future publica- 
tion [33l ]. 
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